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Abstract We investigate masses of light mesons from a coupled system of Dyson-Schwinger (DSE) and Bethe- 
Salpeter equations (BSE), taking into account dominant non-Abelian, sub-leading Abelian, and dominant pion 
cloud contributions to the dressed quark-gluon vertex. The axial-vector Ward-Takahashi identity preserving 
Bethe-Salpeter kernel is constructed and the spectrum of hght mesons calculated. Our model goes significantly 
beyond the rainbow-ladder. We find that sub-leading Abelian corrections are further dynamically suppressed, 
and that our results supersede early qualitative predictions from simple truncation schemes. 
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1 Introduction 

Using the Dyson-Schwinger and Bethe-Salpeter 
equations, there have been successful descriptions 
of pseudoscalar and vector mesons '^•—'—•—•^ using 
the simple rainbow-ladder (RL) truncation. This 
is not unexpected since for the spectrum of light 
pseudoscalar mesons, it is the axial- vector Ward- 
Takahashi identity that governs the pattern of dy- 
namical symmetry breaking. These RL investiga- 
tions have also been recently extended to describe 
baryons t^'^; a status report is given in these pro- 
ceedings ^. 

Considerable effort to go beyond RL, using the 
prescription of l^'iSJ^ bas been made in a number of 
works tifl' Ji'iS'i^' JJ . These however considered sub- 
leading Abelian corrections rather than the dominant 
non-Abelian parts of the quark-gluon vertex ^^'-^ . In 
Ref. ^ we performed a first calculation of the lead- 
ing non-Abelian diagrams in a coupled DSE/BSE ap- 
proach. Here we include the sub-leading Abelian cor- 
rections together with the pion back-reaction tii'iS. 

2 Dyson-Schwinger Equations 

To solve the DSEs for the quark propagator and 
quark-gluon vertex, shown in Fig. [TJa) and Fig. [IJb) 
we must introduce a truncation scheme. The BSE 
must be consistently truncated, being careful to pre- 



serve chiral symmetry. This is expressed via the axial- 
vector Ward-Takahashi identity (axWTI) and ensures 
that pious are the pseudo-Goldstone bosons. 

The simplest truncation that satisfies this crite- 
rion is that of RL whereby the full quark-gluon vertex 
is replaced by a bare vertex The axWTI preserv- 
ing kernel in the BSE then corresponds to a single 
gluon exchange, re-summed to all orders thus provid- 
ing the 'ladder'. These RL models effectively subsume 
additional vertex corrections from the Yang-Mills sec- 
tor and unquenching via an effective coupling. 

To make an extension to the rainbow-ladder trun- 
cation, we investigate the quark-gluon vertex and de- 
termine the impact of corrections beyond tree-level to 
our quarks and mesons UL- 21' ^. ^ , Follow- 
ing the analysis of ti^'^ we approximate the full DSE 
Fig. [IJa) with the (nonperturbative) one-loop struc- 
ture of Fig. [Ijb). Here the first 'non-Abelian' loop- 
diagram in Fig. [ijb) subsumes the first two diagrams 
in the full DSE to first order in a skeleton expansion 
of the four-point functions. We neglect the two-loop 
diagram in the full DSE Fig. [Ija), which is justified 
for small and large momenta ti^'J^. The remaining 
'Abelian' contributions are split into the non-resonant 
second loop-diagram in Fig.[ljb) and a third diagram 
containing effects due to hadron back-reactions. The 
details of including the pion-backreaction are con- 
tained within Refs. see therein for details. 
This provides an extension of our previous truncation 
involving only the non-Abelian diagram ^ . 
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Fig. 1. DSEs for: (a) fully dressed quark propagator; (b) full quark-gluon vertex; (c) truncated quark- 
gluon vertex. Internal propagators are dressed, with gluons shown by wiggly lines, quarks by straight lines 
and dashed lines mesons. White-filled circles show meson amplitudes whilst black-filled represent vertex 
dressings. 



3 Bethe-Salpeter equation 

The Bethe-Salpeter equation describing a rela- 
tivistic bound-state of mass M is calculated through 

Jk 

(1) 

Here = / ^gr, T{p;P) = T^''\p;P) is the Bethe- 
Salpeter vertex function of a quark-antiquark bound 
state. It is a homogeneous equation, with a discrete 
spectrum of solutions at momenta P^ 



-M? 



responding to A(P^) = 1. The lightest of these Mi 
pertains to the ground state solution. The momenta 
fc+ = k+rjP and fc_ = k~{l—r])P are such that the total 
momentum P of the meson is given by P = k^—k^ and 
the relative momentum k = (k^ + k-) /2. The mo- 
mentum partitioning parameter is 77, of which physi- 



cal observables are independent. For convenience we 
choose it to be 1/2 without loss of generalisation. 
The object Kl^{p,k;P) is the Bethe-Salpeter kernel, 
whose Latin indices refer to colour, flavour and Dirac 
structure. 

The Bethe-Salpeter vertex function r'^^(p;P) can 
be decomposed into eight Lorentz and Dirac struc- 
tures. The structure is constrained by the transfor- 
mation properties under CPT of the meson ^ . 

The consequences of chiral symmetry can be 
expressed through the axial-vector Ward-Takahashi 
identity for flavour non-singlet mesons 

-tP^Tl^S^\p+hs+l5Sp\p-)-2mnT\p;P), (2) 

with p and P the relative and total momentum of the 
meson respectively, p± is combination p±P/2, F^^ the 
axial-vector and F^ the pseudoscalar vertex. 
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Fig. 2. The fiavour non-singlet axWTI written showing the connection between the quark self-energy and 
Bethe-Salpeter kernel. All propagators are dressed, with wiggly and straight lines for gluons and quarks 
respectively. The crossed circle indicates a 75 insertion. 
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It is well-known that one may construct a Bethe- 
Salpeter kernel K^^ satisfying the axWTI by means of 
a functional derivative of the quark self-energy. Ap- 
plying this cutting procedure to the quark DSE of 
Fig. [Ija) with the quark-gluon vertex of Fig. [Ijc) 
yields the Bethe-Salpeter equation portrayed in 
Fig. El 12^. Because we do not trivialize the momen- 



tum dependence of our propagators, the calculation 
is genuinely two-loop. Such a truncation allows us 
to consider the dominant non-Abelian corrections ^ , 
together with the sub-leading Abelian corrections ^ . 
We are now in a position to perform a systematic 
comparison of the impact of each vertex correction 
on the spectrum of light mesons. 
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Fig. 3. The axWTI preserving BSE corresponding to our vertex truncation. All propagators are dressed, with 
wiggly and straight lines showing gluons and quarks respectively. 



4 Results 



Model 
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RL 


138 


645 


758 


926 


912 


NA 


142 


884 


881 


1056 


973 


AB 


137 


602 


734 


889 


915 


AB-fNA 


142 


883 


878 


1052 


972 


NA-fPI 


138 


820 


805 


1040 


941 


PDG '^"^l 


138 


400-1200 


776 


1230 


1230 



Table 1. Masses for a variety of mesons cal- 
culated using rainbow- ladder (RL), additional 
corrections from the Abelian (AB) and non- 
Abelian (NA) diagrams, and with the pion 
backreaction (PI). Masses are given in MeV. 

We present results in Table [TJ The benchmark 
for comparison, for which the parameters of the in- 
teraction are fixed to meson observables, are those of 
the rainbow-ladder (RL) approximation. By includ- 
ing the dominant non-Abelian vertex correction (NA) 
we find only a small change in the pion mass, as ex- 
pected since it is a pseudo-Goldstone boson. For the 
remaining bound states we typically see an enhance- 
ment of the mass of the order of 100-200 MeV, de- 
pending upon the meson channel. That is, inclusion 
of the non-Abelian diagram is repulsive in the meson 
channels considered here. 

If we turn off the dominant non-Abelian correction 
and instead consider the sub-leading Abelian diagram 



(AB) we obtain the results labelled AB in Table [TJ 
Based on the relative strength of the correction due 
to its colour factors we expect the results to be 
suppressed, and indeed this is the case. The mass 
of the pion is protected by chiral symmetry and so 
it receives negligible contributions from such correc- 
tions beyond rainbow-ladder. For the scalar, vector 
and Oi axial- vector we see mass reductions of the or- 
der of 20 — 40 MeV, while for the bi there is a slight 
repulsion of 3 MeV. This gives a strong indication 
that such Abelian corrections to rainbow-ladder are 
generally small and attractive. 

Now, if we consider both the Abelian and non- 
Abelian corrections to the quark-gluon vertex to- 
gether we might expect the effects on the meson 
bound-state to stack. That is, for the vector meson 
we would expect to see a ^ 120 repulsion from the 
NA correction, with a ~ 25 attraction from the AB 
correction, resulting in a bound-state mass of ~ 860. 
Instead, we see that the Abelian corrections are heav- 
ily suppressed by the non-Abelian ones for all meson 
channels, giving results that are almost identical to 
the those from the NA diagram alone. The conclusion 
is that naively adding and subtracting the results of 
different independent studies without performing the 
combined dynamical calculation can be misleading. 
It is thus reasonable, with the interaction model pre- 
sented here, to ignore the Abelian diagram completely 
whenever we take the dominant non-Abelian diagram 
into account. 
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Finally, having demonstrated that the Abelian di- 
agram is significantly suppressed in the meson spec- 
trum, we consider the dominant non- Abelian correc- 
tions beyond rainbow-ladder with unquenching effects 
in the form of pion exchange (PI). This gives rise to 
the fifth row of Table [TJ As expected, the inclusion 
of a pion exchange kernel is generally attractive, with 
an 80 MeV reduction of the vector meson mass with 
respect to the RL-I-NA result. This demonstrates, as 
suspected in tiS'iS the near cancellation of beyond- 
the-rainbow corrections with unquenching effects in 
this channel. Here, this is not an exact mechanism 
but the result of dynamical combinations of attrac- 
tive and repulsive components of the quark-gluon ver- 
tex; a more accurate picture will be revealed when 
improved approximation schemes are employed that 
permit quantitative study. 

5 Discussion and Outlook 

We improved upon previous beyond-the-rainbow 
investigations, which focused upon pion unquenching 



corrections til'— and the dominant non- Abelian 
corrections to the quark-gluon vertex ^ by includ- 
ing the sub-leading Abelian corrections to the quark- 
gluon vertex. 

The study we presented involves many techni- 
cal developments in the calculation of the quark 
propagator and quark-gluon vertex at complex mo- 
menta. For highly non-trivial two-loop diagrams such 
as the crossed-ladder kernel, new and improved solu- 
tion methods were required to render the calculation 
tractable. Despite these developments, the study re- 
mains qualitative and serves only as an exploratory 
study of the relevance of the different contributions 
beyond rainbow-ladder considered thus far. 

Future investigations will focus upon including re- 
alistic input from the ghost and gluon sector of the 
theory. These are currently under investigation and 
the results will be reported elsewhere. 

This work was supported by the Helmholtz- 
University Young Investigator Grant No. VH-NG- 
332 and by the Helmholtz International Genter for 
FAIR within the LOEWE program of the State of 
Hesse. 
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